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1 .  INTRODUCTION 


This  technical  report  is  a  brief,  but  significant,  extension  of  an  earlier  report 
Rachele  and  Tunick  (1992)  in  which  the  focus  of  the  study  was  to  determine  the 

sensitivity  of  the  optical  parameter  Cj. j  due  to  variations  of  the  observed 

atmosphere  from  a  true  (ideal)  average  state.  Estimates  of  are  strongly 
dependent  on  the  vertical  gradient  of  the  real  index  of  refraction,  which,  in 
turn,  is  very  sensitive  to  the  vertical  gradients  of  temperature  and  moisture. 
Variations  of  the  atmosphere  from  the  true  average  state,  which  we  characterize 
by  Monin-Obukhov  similarity  equations,  were  envisioned  to  be  due  to  sensor 
response  errors,  sensor  misplacement,  and  data  processing  (resulting  in  incorrect 
values  of  atmospheric  temperature,  pressure,  relative  humidity,  and  wind)  and 
natural  variations  of  the  atmosphere  due  to  naturally  uneven  terrain  and 
nonsteady  state  and  nonhomogeneous  conditions.  We  simplified  our  concept  by 
assuming  that  all  natural  variations  and  errors  due  to  sensor  performance  and 
positioning  could  be  lumped,  and  that  collectively  they  are  normally  distributed, 
even  though  the  natural  variations  are  known  not  to  be  normal  (Pries  and  Appleby, 
1967;  Izumi  1971). 

The  purpose  of  this  study  was  to  approximate  variations  in  the  receiver  coherence 

diameters  and  modulation  transfer  function  (MTF)  due  to  variations  in  .  Many 
of  the  formulations  and  results  presented  in  this  report  were  also  used  in  the 
earlier  report  (Rachele  and  Tunick,  1992). 

We  alert  the  reader  to  a  possible  pitfall  of  misinterpretation  and  application. 
The  so-called  natural  variations  used  in  this  study  were  not  determined  from 
field  data,  instead  they  are  creations  on  our  part  of  what  we  felt  were 
reasonable.  In  particular,  we  not  only  required  that  the  distributions  of  the 
fluctuations  be  normally  distributed,  but  we  also  specified  the  values  of  the 

variances.  (As  such,  the  results  are  at  best  representations  of  ,  rD,  and  MTF 
to  the  input  parameters,  as  specified,  and  should  not  be  interpreted  as  real- 
world  results.)  The  mean  values,  however,  were  determined  from  field  data 
(Stenmark  and  Drury,  1970). 

Critical  to  this  study  is  the  understanding  that  we  deal  with  long-term  (30  min) 
mean  values  and  departures  from  the  long-term  means.  Ideally,  we  assume  that  the 
atmosphere  is  in  steady  state  and  horizontally  homogeneous.  We  recognize, 
however,  for  various  reasons  it  is  not.  We  attempt  to  observe  (measure)  and 
model  the  atmosphere  in  this  ideal  framework,  but  allow  for  adjustments  by 
randomizing  the  basic  inputs  to  the  model.  If  indeed  the  ideal  were  possible, 
then  the  long-  and  short-term  mean  values  that  we  use  would  be  precise  and 
representative  of  the  atmosphere.  We  know,  however,  that  this  picture  is 
unrealistic  to  one  degree  or  another.  Our  concern  then  is  to  determine  the 
combined  effect  of  the  natural  variations  and  sensor  errors  on  the  character  of 
the  atmosphere  and  subsequently  effects  on  the  parameters  of  interest. 

The  types  and  numbers  of  atmospheric  measurements  needed  to  characterize  the 
atmosphere  vary  with  the  models  used  and  with  the  researchers  preferences  and 
budget.  In  this  study,  we  assumed  that  measurements  of  temperature,  pressure, 
windspeed,  and  relative  humidity  are  made  at  two  levels  in  the  lower  atmosphere. 
Each  of  these  basic  measurements  is  averaged  over  a  finite  time  period,  30  min. 
We  will  not,  in  this  report,  discuss  the  sensors,  their  response  characteristics, 
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etc.,  or  the  methods  for  recording  or  processing  the  data.  Instead,  we  assume 
that  the  long-term  (30  min)  values,  however  determined,  can  be  in  "error."  We 
also  assume  that  the  sensors  may  not  be  precisely  positioned  with  height.  In 
what  follows,  we  present  a  methodology,  including  formulations,  for  estimating 

the  sensitivity  of  ,  t0  ,  and  MTF  due  to  the  natural  variations  and  to  sensor 
errors . 

2 .  EQUATIONS 


In  this  section,  we  present  the  basic  expression  for  based  on  Tatarski 
(1961),  which  is  strongly  dependent  on  the  vertical  gradient  of  the  real  index 


of  refraction 


.  Moreover,  has  been  shown  to  be 

dz  dz 


a  function  of  the 


and  specific  humidity, 
dz  QZ 


vertical  gradient  of  potential  temperature, 

If  one  assumes,  as  we  do,  that  the  atmosphere  is  in  steady  state  and  horizontally 

ci0  dci  cLn 

homogeneous,  then  and  and  consequently  can  be  replaced  by  their 


partial  derivatives, 


30  dq 


dn 


30 


,  and  .  In  our  approach  we  represent  and 


and  the  wind  gradient  by  "similarity"  expressions  given  in  equations  (2) 
through  (4)  .  Observation  of  these  equations  shows  that  these  gradients  are 

primarily  functions  of  height  Z,  U* ,  0*,  q*  and  L,  and  hence  the  gradients  are 
a  function  of  the  variations  of  these  quantities.  We  then  proceed  to  show  how 

one  can  formulate  variations  in  u* ,  0*,  q* ,  and  L,  that  is,  du* ,  d0*,  dq*  , 
and  dL.  as  a  function  of  variations  in  V,  P,  T,  Z,  and  f  at  two  heights. 


Next  we  show  that  the  basic  form  of  Cf,  can  be  recast  in  terms  of  0*  and  q* , 
equation  (65).  Hence,  if  0*  and  q*  are  allowed  to  vary  in  a  random  way,  we  can 
compute  variations  of  . 


The  basic  formulation  we  use  for  computing  Cn  is  based  on  Tatarski  (1961)  ;  that 
is, 


C*  =  2> 


(1) 


where 


b  ~  a  constant  -  3.2  (Wyngaard,  1973;  Hill,  1989;  Andreas,  1988) 

u  *  kz 

Kh  -  — t —  =  turbulent  exchange  coefficient  for  heat 
Vh 
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u*  -  friction  velocity 
k  -  von  Karmen's  constant  (0.4) 
z  -  height  above  ground 


,  /  z\1/2 

<pH  =  II  -  15—  j  =  dimensionless  lapse  rate  for  heat  in 


unstable  conditions. 


(\  •* 

4>w  -  — J -  energy  dissipation  rate,  Panofsky  (1968) 

/ 

=  I1  “  15  —  J  =  dimensionless  lapse  rate  for  momentum 

L  -  Obukhov  (1946)  length 
n  -  real  index  of  refraction 


-  height  derivative  of  n 


In  addition,  ——  is  written  as  a  function  of  the  height  derivatives  of  potential 

dz 

cf0  dcr 

temperature  and  specific  humidity,  that  is,  -r-  and  (Tunick  and  Rachele, 

az  dz 

1991).  One  way  of  evaluating  and  is  to  assume  that  they  are  express- 

dz  dz 

ible  in  similarity  form,  which  implies  that  we  are  considering  mean  values.  In 
particular,  for  unstable  atmospheric  conditions,  we  write 


s  -  ■£  (*  -  “ir 


Is .  -2i/i  -  i54r1/2 


dz  kz\ 


z\-V* 


where 


q  -  specific  humidity 
v  -  windspeed 
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0*  -  temperature  scaling  length 

q*  -  specific  humidity  scaling  length 

u*  -  friction  velocity 


L  = 


U*2  T 

*•*  VI 

kg  0t 


(5) 


where 


g  -  acceleration  due  to  gravity 

Tr  -  temperature  at  a  reference  height  Zr 

(1  +  0 . 61qr)  (6) 

(7) 


dO  .  dq  , 

— —  and  — “  ,  we  have 

dz  dz 

imposed  the  conditions  of  steady  state  and  horizontal  homogeneity  on  the 
atmosphere.  Furthermore,  the  similarity  expression  only  makes  sense  in  terms  of 
mean  values,  which,  for  unstable  conditions,  involves  time  averages  on  the  order 
of  30  min. 

As  mentioned  earlier,  the  minimum  amount  of  data  we  require  to  evaluate  equations 
(2)  through  (4)  are  measurements  of  temperature,  pressure,  windspeed,  and 
relative  humidity  at  two  heights.  These  data  are  then  averaged  over  the  total 
time  interval. 

Also,  as  stated  previouslv,  we  assume  that  there  are  errors  in  the  fidelity  of 
measurement,  the  placement  of  the  sensors,  and  the  techniques  for  processing  the 
data,  resulting  in  errors  in  the  estimates  of  the  mean  values.  These  errors  are 
compounded  by  the  natural  transients  that  amplify  the  sensor  errors.  To  quantify 
the  effects  of  these  errors  (variations)  on  CjJ  ,  we  first  integrate  equations 
(2)  through  (4),  resulting  in 


e2  =  8, 


<?2  =  <2i 


ln 

1  y» 

(8) 

+  -21  - 

Lf/7 " x) 

(9) 

k 

l  \y' +  iJ 

v* 

Tvr  -  virtual  temperature  at  reference  height  -  Tr 


0*  =  0*  +  0 . 610rg*  . 


Note  that  to  use  equations  (2)  and  (3)  to  approximate 


8 


( 15c) 


9 


*  1  . 


( 15d) 


Furthermore,  since 


■  (‘  -  ■"!)' 


4*2  *  sz  dz2  *  -rrdL 


where 


■  3  t1  - 


(18a) 


d*a  _  15  f.  _  -.cjjll 

9t  “  4L!l  tj 


(18b) 


Similarly  for  x  -  xt 


where 


,  dx,  ,  dx,  , 

**  =  a^dZl  +  “5^  ' 


a*i  _  _  15f.  _  lsiiV 
dz.  ~  4A1  15  L  ) 


(20a) 


•  J^(l  -  15^1)-^ 
4  L2  L 


(20b) 
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Substituting  equations  (15)  through  (20)  into  equation  (14)  gives 

dv2  -  dv.  *  (Vg  ~  du*  + 

2  1  u* 

15U*  f  _  zi 

kL2  |x2(Xa  -  1)  -  1) 

-  15u*  J  dz2  _  dzi  | 

L  -  1)  ^(X?  -  1)J 


(21) 


(23b) 


(23C) 


_  ^a  ~  (23d) 

30*  ~  0* 
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(24a) 


Since  y  = 


■I1' 


15^/* 

L  j 


where 


Hence , 


and 


dy  _  _  15/,  _  1s-_r\"l/2 
dz  '  2L\  3  L) 


15 

2yL 


dy  _  15  z 

dL  2  yL2 


(24b) 


dy2 


15 

2  y2L 


dz2 


+ 


15  Z2 
2y2L2 


dL  , 


(24c) 


15 

2yxL 


dz2 


+ 


15zt 

2yxL2 


dL  . 


(24d) 


From  equations  (22)  through  (24),  we  obtain 


de-  - 

2  1  e*  kL2 


158*  g2 

kL*  y2  <y|  - 


i)  yx(>i  -  i) 


158* 

kL 


dZr, 


dz2 


(25) 


y2(y;  -  i)  yx(yr  -  D 


To  evaluate  equation  (2b),  we  require  values  of  d8x  and  d82  .  These  expressions 
are  approximated  by  using  the  potential  temperature  expression  for  damp  air, 


8  = 


k'  ~  0.286 


(26) 
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The  differential  of  0  using  equation  (26)  is 


d0  =  0-^  +  0  k'4§  . 

T  P 


(27) 


Therefore, 


*.  •  *  *■* ) 


(28a) 


and 


*■  •  •  **) 


(28b) 


From  equation  (5) ,  we  write 


dL  =  -|*Ldu*  +  -jj=-dTv  +  J^db*v  , 
5u*  dTVx  v*  ae; 


(29) 


where 


dL  m  2  L 
du *  U* 


(30a) 


dL  _  _L_ 

^  "  rVi 


(30b) 


dL  _  _L 

ddi  0; 


(30c) 
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However,  since 


TVi  =  T^l  +  0.61^) 


dT„  = 


ar„ 


dTi 


’idT^ 


ar„ 


d<h 


-dqx 


where 


dTi 


(1  +  0.61^) 


dT 

•*^1 

agx 


0 . 6 1  Tj  . 


To  evaluate  equations  (31)  and  (32),  we  require  values  for  qx  and  dqx. 
obtained  by  using  the  approximations  (Rogers,  1979) 


q  = 


3.8 


f  exp 


l 

273.16 


dq 


Pdp  +  +  4jdr  , 

op  df  dT 


where 


iff  =  -2 

dp  P 


iff  =  st 

df  f 


iff  ,  sss 

dT  T8 


(31) 

(32) 

(33a) 

(33b) 

These  are 

(34) 

(35) 

(36a) 

(36b) 

(36c) 
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f  -  relative  humidity  in  decimal  form 
K  -  5.44  x  103  in  c . g . s .  units 

Substituting  equation  (36)  into  equation  (35)  for  z  -  zx  gives 


d<?i  -  -4t4pi  +  Tdfl  + 

*i  ri  if 


where  qx  is  evaluated  by  using  equation  (34). 

To  evaluate  equation  (32),  we  also  require  values  of  0^  and  dO^. 
obtained  using  expressions  presented  in  Rachele  and  Tunick  (1991) 

K  =  0*  +  0.611^*  , 


and 


dO^  =  d0*  +  O.eiTjdg*  , 


where  q*  and  dq*  are  discussed  in  section  2. 

Substituting  equations  (30)  through  (39)  into  equation  (29)  gives 

dL  =  ^  du* 
u* 


+ 


+  0.61gx)  dTx  +  0 . 6 1  Tx 


+ 


(d0*  +  0.6ir,dg*)  . 

e: 


From  equation  (9),  we  write 


dq  -  l£d,‘ * 


(37) 

These  are 

(38) 

(39) 

(40) 

(41) 
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where 


However, 


where 


(42a) 


_dg_  a  g>  ~  (42b) 

a<7*  <7* 


(42c' 


(42d) 


since 


y' 


-  15 


z\l/i 

LI 


(43) 


dy'  -  &-dz  ♦  &dL  , 


(44) 


SxL  «  —12- 

dz  2Ly' 


(45a) 


dy '  _  15z 

dL  2h\y' 


(45b) 
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Substituting  equations  (42)  through  (45)  into  equation  (41)  gives 


dg2  -  dqx  = 


-«±dq-  ♦  ISfllf _ - - 5*8 - 1 

r*  *L  | A(y?  - 1)  yafya  -  i>  J 


+ 

kL 2 


- h - - h - \dL  , 

1  |ya(y22  -  i>  yi(yi2  -  D  J 


where 


(dp,  df,  a.fc ,  dp,  df,  Jcdr,  I 

-^-pr  *  +  if  dT»  *  ‘  '  ai_jTj  ' 


At  this  point,  we  have  four  equations,  (21),  (25),  (40),  and  (46),  in  four 
unknowns  that  are  driven  by  deviations  in  V,  T,  P,  z,  and  f  shown  in  matrix  form. 
This  is  our  solution  matrix. 


Yii  0  0  Yi4  du* 

0  Y22  0  Ya*  d0* 

0  0  Y33  Y34  dg* 

Y«i  Y42  Y43  Y44.  .  ^ 


where 


( v2  -  Vx) 


iU*f  z2 _ Zl 

L2  jx^(x£  -  1)  xx(xj  -  1) 
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Ax  =  dv2  -  dvx  + 


15  u' 


dzx 


Y24  = 


Y34  - 


“V1  ' 

***  \*a(*2  ~  1)  *x(xi  - 

_  150* 

za  Zl 

kL2 

y2(ya  -  1)  yx(yl  -  1) 

„  .  <e2  -  0x) 

-  e. 

■+? 

■  *'*)  -  Kf  -  *'x) 

150*  f 

dz2  dzx  \ 

kL  ly2(y22  -  i)  y^yi-i)/ 

„  -  <?2  -  <*1> 

T“ ■  «• 

15g* 

Z2  _  Z1 

icL2  ' 

y2  (y22  -  1)  yi(yia  -  1) 

(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


"a 

*2 


(57) 


+ 


1 5  cj*  dz2 

kL  U(y?  - 1) 


dzx 

y[  (y!!  -  1) 
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(58) 


2L 

Y“  =  ^ 


L  L 

Y“  * '  <F  '  Tt 


Y43  *  -0.61 - ^ 

43  0: 


Y44  “  "I 


A4  *  -LdTVi 


dTVi  *  (1  +  0.6^)  +  0.612^4?! 


dq±  =  ~q: 


(dp1  dfx  K 

\ p.  a  if 


drj  . 


The  matrix  equation  (48)  is  solved  for  randomized  sets  of  deviations  of  T,  P,  V, 
f,  and  sensor  positions,  resulting  in  distributions  of  du*,  dL,  d9*,  and  dq*  . 

Simultaneously,  we  compute  for  each  sample  set  resulting  in  a  distribution 

of  .  The  formulation  we  use  for  evaluating  Cjj  (see  Rachele  and  Tunick,  1991) 
is 


where 


c*  =  A'e*a  +  fl'ev  +  Cq *2  f 


A'  =  £(6.241  X  10"9)  j£lT2/3Z-2/3 |l  -  Y-|j  ^  ) 


19 


(67) 


B1  =  jb  [3 . 11  x  10'9]  -~lr2/3z'2/3|l  -  y-5p{  } 


C'  =  i?(3 . 88  x  io-10)  -~ir2/3z-2/3|i  -  ) 


(68) 


{}  = 


'69) 


where 

b  -  3.2 
P  =  Y  =  15 

3.  RECEIVER  COHERENCE  DIAMETER  (ro) 

The  formulation  we  use  for  computing  the  receiver  coherence  diameter  is  discussed 
in  Miller  and  Ricklin  (1990), 

r0(D  =  |o.423ic2^C^(s)  ^(S)  dsj  3/5  ,  (70) 


where 


ds  -  a  path  length  increment 
r  -  propagation  path 
K  -  wave  number 

/  S\5/3 

Qj(s)  -  (  —  1  for  a  spherical  wave 

L  -  the  total  path  length 
4.  THE  MODULATION  TRANSFER  FUNCTION 

The  modulation  transfer  function  (MTF)  selected  for  this  study  is  also  discussed 
in  Miller  and  Ricklin  (1990).  In  particular,  Miller  and  Ricklin  refer  to  it  as 
the  near-  and  far-field  slow  turbulence  modulation  transfer  functions  (TMTF) , 
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TMTF1(f,  r) 


>|-3 . 44  ~ 


where  R  is  the  focal  length,  A  is  the  wavelength,  and  f  is  the  spatial 
frequency . 

5.  EXAMPLES 

For  this  study,  we  selected  two  sets  of  data  from  Davis,  California,  during 
unstable  conditions  with  clear  skies  (Stenmark  and  Drury,  1970).  The  first  case 
was  chosen  to  emphasize  the  effects  on  atmospheric  stability  of  moderate  (5.0  to 
7.0  m/s)  winds.  The  second  case  contains  lighter  (2.0  to  3.0  m/s)  winds  and 
somewhat  intensified  atmospheric  instability.  The  Davis  field  site  (a  flat,  5-ha 
area  at  17  m  above  sea  level)  is  located  about  2  km  west  of  the  main  portion  of 
the  University  of  California  at  the  Davis  Campus,  24  km  west  of  Sacramento  and 
113  km  northeast  of  San  Francisco.  The  data  were  taken  during  periods  when  the 
surrounding  fields,  for  the  most  part,  were  crop  covered  and  well  irrigated, 
giving,  in  effect,  homogeneous  surface  conditions  with  respect  to  temperature  and 
moisture.  Advection  effects  were  considered  to  be  negligible.  Profiles  of  wind, 
temperature,  and  specific  humidity  were  measured  at  nine  levels  from  25  to  600 
cm.  Raw  data  were  processed  to  give  one -half  hour  average  profiles.  Table  1 
gives  the  two -level  (1  and  10  m)  values  for  windspeed,  temperature,  specific 
humidity,  relative  humidity,  and  pressure  for  the  two  data  sets  discussed  above. 
From  these  data  we  found  values  of  the  friction  velocity  u*.  the  potential 
temperature  and  specific  humidity  scaling  parameters  8*  and  q*.  and  the  Obukhov 
lengths  L  and  L*,  which  are  also  given  in  table  1. 

TABLE  1.  MICR0METE0R0L0GICAL  PARAMETERS  FROM  THE  DAVIS,  CA,  DATA 


6 .  RESULTS 


As  in  our  earlier  report,  we  caution  the  reader  to  a  possible  pitfall  of 
misinterpretation  and  application  of  the  data  presented  in  this  report.  The  so- 
called  natural  variations  used  in  this  study  were  not  determined  from  field  data; 
instead  they  are  creations  on  our  part  of  what  we  felt  were  reasonable.  In 
particular  we  not  only  required  that  the  distributions  of  the  fluctuations  be 
normally  distributed,  but  we  also  specified  the  values  of  the  variances.  (Note 
however  that  the  mean  values  were  determined  from  field  data- -see  section  5.) 
The  standard  deviation  for  the  windspeed  distributions  was  specified  to  be  3.33 
cm/s.  We  did  not  try  to  adjust  the  windspeed  variance  to  changes  in  the  magni¬ 
tude  of  the  windspeed  itself.  The  standard  deviation  of  temperature  was  chosen 
to  be  0.033  K.  The  standard  deviation  of  relative  humidity  was  one -half  of  1 
percent  relative  humidity,  and  the  standard  deviation  for  the  pressure  was  0.33 
mbar . 

For  each  case  we  computed  distributions  of  U* ,  0*,  q* ,  and  L.  We  then  computed 
distributions  of  (fn  and  rD  (receiver  coherence  diameter).  Finally,  we  derived 
the  near-  and  far-field  slow  MTF  from  the  minimum,  maximum,  and  mean  rQ  and 
plotted  the  slant  path  case  (1  km  length  by  20  m  height)  for  a  beam  pointed 
upward  to  a  target  and  for  a  beam  pointed  downward. 

The  results  of  the  two  cases  presented  in  this  section  are  at  best  representa¬ 
tions  of  the  sensitivity  of  and  subsequently  r„  and  MTF  to  the  input 
parameters,  as  specified,  and  should  not  be  interpreted  as  real  world  solutions. 

Figure  1*  shows  the  two-level  model  (Case  1)  normal  distributions  for  H,  L' E,  U*, 
and  Vr  and  is  presented  to  show  the  fidelity  of  the  normal  distributions 
generated  and  used  as  input  parameters. 

Figure  2  shows  the  derived  distributions  for  T*  (0*)  ,  Q* ,  and  L.  The  range  for 
these  distributions  is  approximately  a  factor  of  2.  Note  that  for  this  case  one, 
since  U*  is  relatively  large  and  T*  is  small,  the  mean  value  for  the  Obukhov 
scaling  length  is  quite  large,  indicative  of  weakly  unstable  atmospheric 
conditions . 

Figure  3  illustra  es  the  random  distribution  for  resulting  from  variations 

in  u* ,  0*,  q* ,  L,  and  Vr.  The  range  for  the  random  distributions  of  is 

approximately  a  factor  of  3.  Additionally,  a  cumulative  distribution  for  is 

shown.  It  suggests  that  about  47  percent  of  the  time  will  have  a  value  equal 
to  or  less  than  its  mean  value. 

Figures  4  and  5  show  the  distribution  of  r0  and  its  effect  on  the  MTF  for  Case  1. 
Note  that  the  variations  of  the  MTF  are  based  on  extreme  values  of  r0.  Figure 
4  gives  the  MTF  for  a  slant  path  looking  upward,  and  figure  5  shows  the  MTF  for 
a  slant  path  looking  downward. 


♦Figures  are  presented  at  the  end  of  the  text. 
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Figures  6  and  7  show  the  distributions  of  Vz,  U* ,  T* ,  q*  ,  and  L  for  Case  2,  and 

figure  8  shows  C%  for  Case  2.  Note  here  that  mean  values  for  T*  and  q*  are 
slightly  greater  in  magnitude  than  those  from  Case  1,  and  the  mean  value  for  U* 
is  small.  This  results  in  L  values  smaller  in  magnitude  for  this  case, 

representing  more  moderately  unstable  atmospheric  conditions.  is  signifi- 
cantly  larger  in  magnitude  (that  is,  on  the  order  of  10"13)  and  its  range  is 
approximately  a  factor  of  5. 

7.  SUMMARY  AND  CONCLUSIONS 

The  purpose  of  this  study  was  to  examine  the  sensitivity  of  the  receiver 

coherence  diameter,  rol  and  the  MTF  for  deviations  of  Cjj  due  to  variations  of 
the  natural  atmosphere  and  from  misrepresentations  of  the  atmosphere  due  to 
errors  in  measurement  and  processing  of  the  basic  parameters  involved.  We 
selected  two  sets  of  data  from  Davis,  California,  taken  during  unstable 
conditions  under  clear  skies.  We  presented  graphs  illustrating  the  derived 
distributions  of  the  similarity  scaling  parameters,  cj| ,  r0,  and  plots  of  the  MTF 
for  two  slant  path  situations. 


We  found  that  for  the  specified  variations  of  Vr,  T,  P,  and  f  that  could  be 

determined  within  a  factor  of  3  for  Case  1  and  a  factor  of  5  for  Case  2. 


The  question  remains  as  to  what  are  acceptable  ranges  or  values  for  .  The 
answer  lies  wholly  with  their  use,  or,  that  is,  it  depends  on  the  application  for 

.  We  found  for  instance  that  r0,  the  receiver  coherence  diameter,  can  vary 
from  1.26  to  3.76  cm  for  Case  1,  causing,  in  turn,  a  significant  effect  upon  the 
near-  and  far-field  slow  MTF.  For  Case  2  r0  varied  from  3  to  6  and  also  resulted 
in  more  skewed  distributions  of  MTF. 

In  conclusion  we  feel  confident  that  our  methodology  was  a  sound  indicator  of 
potential  variations  in  optical  parameters  in  our  models.  Needless  to  say,  a 
more  thorough  study  should  be  performed  using  more  complete  field  data. 
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Figure  l.  (cont.) 
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Figure  2 .  (cont) 
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4.79  6.83  8.87  7.90  8.94  9.97  11.01  12.04  13.08 

Cn**2  (  NT14) 

Mean  Cn**2  -  7.90  x  1014  (  ffl2'*  ) 


(a) 


A 


Cn**2  (10*1} 

Mean  Cn**2  ■  7.90  x  10'14(  m'2,a) 

(b) 

Figure  3.  Case  1  -  Distribution  for  : 

(a)  random  and  (b)  probability. 
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Transfer  Function 


Rnot  (cm) 

Mean-  3.92  cm,  Path-lkm  x  20m  ht,  0.55um 


Figure  4.  Case  1  -  (a)  Random  distribution  for  r0  and  (b)  the 
resultant  MTF  for  a  slant  path  looking  upward. 


Rnot  (cm) 


Mean-  2.24  cm,  Path-lkm  x  20m  ht,  0.55um 

(a) 


THE  NEAR  AND  FAR  REID  SUOW  TVfTF 
(Toed  length  =  10cm,  wavelength  =  055  micronsj 
[Ro-mm  -  159cm,  Ro-max  =  3.95cm,  Ro-mean  =  Z23crrj 

(b) 


Figure  5.  Case  1  *  (a)  Random  distribution  for  r0  and  (b)  the 
resultant  MTF  for  a  slant  path  looking  downward. 
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Figure  6.  Case  2  -  Random  distribution  for:  (a)  u* 
and  (b)  normal  distribution  for  Vr. 


Figure  6.  (cont) 
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Specific  Humidity  Scaling  Parameter 


Mean  ■  -4.954  x  10  '*  (g/g) 


(b) 


Figure  7.  Case  2  -  Random  distribution  for: 
(a)  T*,  (b)  q*,  and  (c)  L. 
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Figure  7.  (cont) 
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Figure  8.  Case  2  -  Distribution  for  C%: 

(a)  random  and  (b)  probability. 
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Figure  9.  Case  1  -  (a)  Random  distribution  for  r0  and  (b)  the 
resultant  MTF  for  a  slant  path  looking  upward. 
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Figure  10.  Case  1  -  (a)  Random  distribution  for  r0  and  (b)  the 
resultant  MTF  for  a  slant  path  looking  downward. 
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